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Abstract 

A p-spin interaction Ashkin- Teller spin glass, with three independent Gaussian prob- 
ability distributions for the exchange interactions, is studied by means of the replica 
method. A simple phase diagram is obtained within the replica-symmetric approxi- 
mation, presenting an instability of the paramagnetic solution at low temperatures. 
The replica-symmetry-breaking procedure is implemented and a rich phase dia- 
gram is obtained; besides the paramagnetic phase, three distinct spin-glass phases 
appear. Three first-order critical frontiers are found and they all meet at a triple 
point; among such lines, two of them present discontinuities in the order parame- 
ters, but no latent heat, whereas the other one exhibits both discontinuities in the 
order parameters and a finite latent heat. 
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1 Introduction 



In recent years much progress has been achieved in the understanding of mag- 
netic disordered systems. Among those, one may single out the spin glasses 
(SGs) [1-3], for which significant advances were obtained, as a result of a 
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large effort dedicated to them. Most of the SG studies were carried for two- 
spin-intcraction models, that were at the early stage, intended to explain the 
physical behavior of some peculiar magnetic compounds. Nowadays, apart 
from this motivation, such SG models were identified to be closely related to 
a large diversity of physical problems, like neural networks, protein folding, 
and optimization problems [2,3]. The mean-field theory of the Ising SG model, 
considered in terms of an infinite-range- interaction model [4], is nowadays ac- 
cepted as satisfactorily understood. The simplest solution, based on a single 
SG order parameter, known as replica-symmetric (RS) solution [4], presented 
serious difficulties, and it was shown to be unstable at low temperatures [5]. 
The correct mean- field solution for the Ising SG was proposed by Parisi [6], 
and it consists in an infinite number of order parameters - a procedure called 
of replica-symmetry breaking (RSB). However, it is not clear up to the mo- 
ment, whether the mean-field solution is appropriate for the description of real 
- short-range-interaction - SG systems [3] . 

The j9-spin-interaction {p > 2) models were initially introduced as merely the- 
oretical problems [7,8]; however, the identification of a close analogy between 
the transitions occurring in p-spin interaction SG models and those obtained 
in mode coupling theories of structural glasses [9,10], motivated many studies 
on such models. A nice feature of infinite-range p-spin-interaction SG models 
is that in the p — > cxo limit the energy levels become independent random 
variables, yielding an exactly solvable model, known as the random-energy 
model [7,8]. Also, for p > 2 the phase transition scenario of these SG models 
is quite different from the one found in the corresponding p — 2 model. Besides 
the usual equilibrium transition temperature, there exists another transition 
temperature presenting a dynamical character. Right below the equilibrium 
transition, a single step in the RSB procedure is sufficient for stabilility [11], 
i.e., the order parameter is properly represented in terms of a one-step Parisi 
function. For models where there is some kind of competition between different 
types of interactions, the one-step RSB order parameter changes dramatically 
the phase diagrams obtained from RS solutions [12,13]. Usually, in the p ^ oo 
limit, the one-step RSB yields the correct solution, revealing new phases which 
are not present in the corresponding RS approach. 

It should be mentioned that, recently, there has been an increasing interest 
in the study of p-spin-interaction SG models, motivated either by attaining a 
better understanding of such models, or by the possible applications in other 
fields of science. The study of infinite-range ]9-spin-interaction SG models, by 
means of the replica method, has produced many interesting new features 
as a consequence of the inclusion of magnetic fields [14,15], ferro- [16-18] 
and antiferromagnetic [19,20] interactions, as well as a competition between 
quadrupolar and SG orderings [13]. Recent dynamical studies also have been 
carried [21,22], leading, in particular, to an analysis of the barriers separating 
metastable states [21]. Besides that, such models have been investigated lately 
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through rigorous approaches [23] and in a quantum version [24]; new apphca- 
tions were explored, e.g., connections to the protein folding problem [25], to 
error correcting codes [26], and to many biological systems [27]. 

The Ashkin- Teller model [28] is one of most studied systems in statistical me- 
chanics, due mainly to the richness of critical phenomena revealed by its phase 
diagrams both in two and three dimensions [29] . This wealth of results steams 
from the competition between the two- and four-spin interactions present in 
the model. Herein we consider an infinite-range Ashkin-Teller-like SG with 
p-spins interactions in order to investigate the effects of an analogous compe- 
tition. We present the full solution in the p ^ oo limit. The resulting phase 
diagram shows three distinct phases, similarly to what has been found in 
other Ashkin- Teller SG models with p = 2 [30-32]. However, the nature of 
the present transition lines are changed, and we get a triple point common to 
the three phases, instead of a multicritical one as found in the previous p — 2 
works [30-32]. In particular, wc find that in the limit of two independent 
Ising-like models (or 4-state clock model) the equilibrium transition occurs at 
a multiphase point. We also show that, for p — oo, some particular cases of 
the present model are equivalent to random-energy models with uncorrelated 
energy levels. 

This paper is organized as follows. In the next section we introduce the model 
and determine the free-energy density functional, obtained through the replica 
method. In section 3 we determine the phase diagram within the RS solution 
and consider the stability of such a solution against Gaussian fluctuations 
in replica space. In section 4 we apply the first stage of Parisi's Ansatz to 
determine the phase diagram. Our main conclusions are drawn in section 4. 
Equivalences with random-energy models, in the limit p — > oo, are shown in 
an appendix. 



2 The Model 

In the present work we will consider a p-spin interaction Ashkin-Teller-like SG 
model, defined by the Hamiltonian 



i<n<...<ip<Ar 

where a and r (= ±1) are Ising spin variables. All interactions are infinite- 
range- like, and similarly to what has been done in previous works [30,31], 
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herein we will be restricted to the case of independent couphngs J-^^^i^ (a 
1, 2, 4), each one following its own Gaussian probability distribution, 



The free-energy density is given by 



exp < 



(2) 



(3) 



where (3 = (ksT) ^, and ( ) represents an average over the disorder. In the 
following we will use the replica method to calculate 



(InZ) — lim 



(Z") - 1 



n 



(4) 



Performing the averages over the random couplings, one gets 



(Z") =Tr exp 



4iVP- 



l<ii<—<ip<N \a=l 



+ 



+ 



47VP-1 



l<ji<---<ip<Ar \o=l 



(5) 



El \ ^ „a a a a a 



l<ii<—<ip<N \a=l 



where a = 1, . . . ,n represents the replica index. As usual, the sums over p 
sites may be reduced to sums over a single site, e.g.. 



I \ ^ „a „a \ \ ^ \ ^ „a „a Jj „a 



l<h<—<ip<N \a=l / a,b=l l<ii<-- <ip<N 

= ^E (j:->'Y+o{n^- 



P- P- (ab) i / 

where J2{ab) denotes a sum over distinct pairs of replicas. One may now intro- 
duce, for each distinct pair of replicas {ab), the order parameters. 
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Ql,ab 
(l2,ab 
rob 



\ ^ „a„b /„a„b\ 



-E 



a^b I ^a^b 



N 



T-T 



T T }, 



\^ -a_.a„b^b I „a^a „b^b\ 



(6.a) 
(6.b) 
(6.C) 



as well as their respective auxiliary fields (Lagrange's multipliers) 7i,ab, 72,ab, 
and through standard identities, e.g., for the pair (qi,7i). 



where. 



(ab) 

L^-Y.li,ab iNq,,ab-Y.<<'i]- 

(ab) V i / 



(10) 



(11) 



Applying a similar procedure for the pairs of matrices (q2, 72) and (r, ^), one 
may write 



00 flOO 



00 J —lOO 



n d'yi,abdqi,abd^2,abdq2,abd^abdrab 



(ab) 



X exp {-Ng^a^i^ab}, {qi,ab}, {72,ab}, {q2,ab}, {Cab}, {^^ab})] (12) 



where 



fi'n — E (7l,ab9l,ab + 72,ab92,ab + Cab^" ab) 
(a,b) 



-TE(AU + ^2Xab + ^: 



(a,b) 



< E [7l,ab<C 
Uab) 



- In Tr exp { ^ 1 7i,„,<a,^ + 72,abrf r^* + 



a' ab) 



(13) 



Substituting eq. (12) into eq. (4), we obtain the free-energy density functional 



Pf = -^n({7l,ab}, {gi,ab}, {72,ab}, {g2,ab}, {Cab}, {^^ab}) (14) 



5 



where gn stands for the global minimum of Qn, taken with respect to the 
variational parameters. 

One must now look for stable solutions of the saddle-point equations. The 
simplest Ansatz consists in considering all replicas in an equal footing, and 
thus assume the RS solution. In the next section we will analyse this solution, 
and it will be shown that the only stable RS solution is the paramagnetic one. 



3 Replica- Symmetric Solution 

The RS solution is obtained by considering 

qi,ab^qi, ?2,a6 = ?2, and Tab^r, V(a6), (12.a) 



as well as 



7i,a6 = 7i, 72,a6 = 72, and Cab = C, V(a6). (12. b) 



Performing standard simplifications through Gaussian identities, we get the 
free-energy density functional as 



/ = f [JM - 1) + ^2(?f - 1) + J!(r' - 1)] (16) 
[7i(?i - 1) + 72(^2 - 1) + C{r - 1)] - ^((ln5(x, y, z))),,y,, , 



where 



'^{x, y, z) = 4[cosh(v^x) cosh(^|/) cosh(y^2;) 

+ sinh(^a;) sinh(^|/) sinh(y^^)]. (17) 



In Eq. (16) the double brackets {{))x,y,z stand for Gaussian averages with 
respect to the set of variables {x, y, z), e.g., for an arbitrary function ip{x, y, z), 
one has 



// / \\\ f°° f°° f°° dxdydz 1.^2 , 2 ,^2^ 

J —00 J —00 J —00 (27r)2 
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The parameters qi, q2, r, 71, 72 and ^ in Eq. (16) may be determined from 
the saddle-point conditions. Thus, the equations relating the auxiliary fields 
(71 ) 72,0 to the order parameters {qi,q2,r) become 

7i = lpiPJifqr\ l2 = \p{PJ2fqt\ i = \p{PJ^fr^-\ (18) 

whereas the order parameters are self-consistently given by 

gi = {{ipl{x, y, z)))^ y^^ , q2 = {{(pl{x, y, z)))^ ^^^ , r = {{cpl{x, y, z)))^^^^^ , 

(19) 

with 



ipi{x,y,z) 
'P2{x,y,z) 
'P3{x,y,z) 



— ^tanh{y^x) + ta.nh{^/J^y) tanh(Y^2;) 
[tanh(^O^y) + tanh(^7r3:^) tanh(Y^2;) 
[tanh(y^z) + tanh(v^a;) tanh{y/^y) 



(17.a) 

:i7.b) 

(17.C) 



and D — 1 + tanh(iy7ix) tanh(iy^y) tax^h{^/^z). 

It would be interesting to investigate the full phase diagram obtained from the 
above equations. However, similar to what has been done in previous studies of 
p — 2 Ashkin- Teller models, the most interesting case is the isotropic one, for 
which J2 — Ji - In such a case, besides the usual spin reversal symmetries, there 
is a new symmetry in which the a and r spins variables may be interchanged. 



For p > 1 there is always a disordered, paramagnetic, phase where qi — q2 — 
r — and 71 = 72 = C = 0? with the free energy given by 



/p = -^(2J| + j|)-2Tln2, (21) 

where T is the temperature (herein we work in units such that k-Q = 1)- From 
the above expression we obtain the entropy density 



sp = -^(2j2 + j|)+21n2. (22) 

For low temperatures this entropy becomes negative and the system gets 
frozen. The paramagnetic phase can exist only above the curve shown in Fig. 
1, which is given by 
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Fig. 1. Phase diagram of the p-spin interaction Ashkin-Teller SG, within the RS 
solution, for arbitrary values of p. At low temperatures, the paramagnetic solution 
(P) becomes unstable, and there are no stable non-trivial solutions. 



T 1 / 72 \ 1/2 

There are many other non-trivial RS solutions. The stability analysis of those 
solutions against rephca fluctuations can be performed along the hues pio- 
neered by de Almeida and Thouless [5] for the p — 2 infinite-range Ising SG 
model. In general, the eigenvalues are roots of cubic equations, and become 
too lengthy to quote their expressions, except in a few particular cases. For 
instance, to the 51 = ?2 = 0, 71 = 72 0, r 7^ solution there corresponds, in 
the limit n — > 0, the following longitudinal 



AL = -(/3J4mp-l)r 



p-2 



1 - -{(3Jifp{p - l)r^'^ (1 - 4r + 2>t) 
2 



(24) 



and transversal 

At = \Wpip - l)r^-' [1 - \w^pij^ - 1)^^"' (1 - 2r + t) 



(25) 



eigenvalues, where 

When r = 0, one gets the paramagnetic solution and, from the vanishing of the 
above eigenvalues, such a solution is marginally stable for p> 2. On the other 
hand, in the p 00 the r = 1 solution presents both eigenvalues negative, 
being completely unstable. The same happens for any other non-trivial RS 
solution. Thus, we must look for RSB solutions; this will be done in the next 
section. 
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4 Replica Symmetry Breaking Solution 



Since there is no stable SG solution in the p — > oo limit, within the RS 
assumption, we need to look for RSB solutions. Fortunately, it is enough to 

consider just the first step in Parisi's RSB procedure in order to get the correct 
solution in such a limit, as verified in other j9-spin-interaction models studied 
before [11-13,33]. In the present case, since we deal with three distincts pairs 
of conjugated matrices, namely (qi, 71), (q2, 72), and (r, ^), one should apply 
the RSB scheme for all six matrices. In principle, one could apply the RSB 
scheme for each matrix in an independent way, i.e., each matrix should have 
its own block sizes. However, it is easy to convince oneselves that this is 
not a physically acceptable procedure, due to the symmetries involving the 
interchanging between the spins in the Hamiltonian. Therefore, herein we will 
apply the same RSB for each matrix, i.e., we will divide all six n x n matrices 
into n/m groups of size ra. Following Parisi [6], we denote the elements of the 
off-diagonal blocks by g'1,05 Q'2,0; ''"O; 7i,0; 72,o and and those of the diagonal 
blocks by 52,1, 71,1, 72,1 and ^1. Thus, we obtain the free-energy density 
functional 



-{Jhli + Jhl. + Jlrl){m - 1) + ^{Jlq% + Jlq% + J^rg) 



1 Tfl 

+ ^(7i,i?i,i + 72,1^2,1 + - 1) - ^(7i,ogi,o + 72,0^2,0 + Coro) 

-^(J^ + J| + J|) + ^(71,1 + 72,1 + 

- — ((In Z{xq, xi, X2))) (27) 



where 

Z{xo, xi, X2) = ((A™(xo, xi, X2, yo, yi, ^2))) 2/0,2/1,2/2 , (28) 

A — A cosh(xio) cosh(iii) cosh(xi2) + 4 sinh(tio) sinh(iii) sinh(ii2), (29) 
with uq, ui and U2 defined by 



uo = ^ffifi + ^71,1-71,0 yo , xi + ^72,1 - 72,0 yi , 

U2 = X2 + - Co y2 ■ (30) 
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By extremizing the above free-energy density one gets the equations of state 
within the one-step RSB procedure (see Appendix A); it is easy to see that 
such equations present several non-trivial solutions in the limit p — > oo. A 
careful analysis of the free energy shows that only four of those solutions are 
physically acceptable. For low temperatures there is always a SGI phase (see 
Fig. 2), in which qi^i = ^2,1 = ri = 1, gi,o = ^2,0 = ro = 0, and 



2V2h[2 T 

m = . — . (31) 

^2+ (74/^2)2 

The corresponding free-energy density is given by 



/, = -^(4J| + 2J|)ln2 , (32) 

which yields a vanishing entropy. The transition between the SGI phase and 
the paramagnetic one occurs along a line given by 



where both phases present zero entropy, i.e., there is no latent heat. This line 
holds as long as J4/J2 < V^- Beyond that point, we find two other first-order 
transition lines confining a SGII phase, as shown in Fig. 2, throughout which 
one has ri = 1, m — 2^/\n2 T/J4, and all other order parameters zero. This 
phase presents a free-energy density 



= -v^ J, _ J. _ Tln2 , (34) 



which leads to the entropy density 



^n{T)^-^,+\n2. (35) 

Due to discontinuities in the order parameters, the transition separating the 
paramagnetic and SGll phases is first-order; this critical frontier may be de- 
termined by demanding the free-energy densities given by Eqs. (21) and (34) 
to be equal, leading to the straight line. 



T = — ^ . (36) 

Along this critical frontier there is also no latent heat. For lower temperatures 
another transition occurs, this time from the SGII to the SGI phase. From 
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2 J /J 4 

4 2 

Fig. 2. Phase diagram of the p-spin interaction Ashkin- Teller SG, within a one-step 
RSB procedure, in the limit p ^ oo. The borders of the paramagnetic phase (P) 
present discontinuites in the order parameters, but no latent heat, whereas the crit- 
ical frontier SGI-SGII is a genuine first-order phase transition, with discontinuities 
in the order parameters and a finite latent heat. The phases P, SGI, and SGII co- 
exist at a triple point (TP), following the standard Gibbs phase rule. A distinct SG 
solution (SGIII) becomes possible at the multiphase point (MP), where the phases 
P, SGI, and SGIII coexist. 

their free-energy densities, given respectively by Eqs. (34) and (32), we find 
the corresponding critical frontier. 



T 
T2 



[2 + V2) 



1 + 



J? 
2 J? 



J. 



' -(l + x/2)- 

2 "^2 



(37) 



It is easy to verify that along the SGI-SGII critical frontier there is a finite 
latent heat. Therefore, in this case one has a genuine first-order transition line. 
The three transition lines given by eqns. (33), (36) and (37) merge together 
at a triple point (TP), located at J4/J2 = V2 and J2/T — y/2\n2, as shown 
in Fig. 2. Although two of these transition lines do not represent conventional 
first-order transitions in the thermodynamic sense, since we do not observe 
any discontinuity in the first derivatives of their corresponding free-energy 
densities, i.e., no latent heat, the picture around the triple point follows the 
standard Gibbs phase rule [34]. 

Besides the solutions discussed above, there are also two equivalent solutions 
(which we call SGIII): qi^i = 1, m = 2\/\n 2 T/J2, with all remaining order 
parameters zero, and its symmetric one, i.e., ^2,1 = 1, m — 2y/\n2 T/J2, 
with all remaining order parameters equal to zero. However, the SGIII phase 
is realized only at the multiphase point (MP), [J4 = 0, T/J2 = l/(2-\/ln2)], 
where it coexists with the paramagnetic and SGI phases. 



11 



5 Conclusion 



We have studied a version of the Ashkin- Teller SG model, with p-spin inter- 
ations, by means of the replica approach. The RS solution leads to a simple 
phase diagram, for arbitrary values of p, with a paramagnetic phase that is 
stable at high temperatures, becoming unstable at low temperatures; within 
such a solution, there arc no stable non-trivial solutions at low temperatures. 
By applying a one-step RSB procedure, we have found a rich phase diagram 
in the p — > oo limit, with four distinct phases, namely, a paramagnetic one at 
high temperatures, and three SG phases (SGI, SGII, and SGIIl) at lower tem- 
peratures. The borders of the paramagnetic phase present discontinuities in 
the order parameters, but no latent heat, whereas the critical frontier separat- 
ing phases SGI and SGll is a genuine first-order phase transition, exhibiting 
both discontinuities in the order parameters, as well as a finite latent heat. 
These critical frontiers all meet at a triple point according to the standart 
Gibbs phase rule [34], similarly to what happens in a previously investigated 
model [13]. The SGlll solution is stable only at a multiphase point, where it 
coexists with the paramagnetic and SGI phases. Also, in the p ^ oo limit, 
it is possible to show the equivalence of the model considered herein with a 
random energy model, as defined by Derrida [7,8], which can be solved by 
other methods; this equivalence is shown in Appendix B. A detailed analysis 
of the corresponding random-energy model will be published elsewhere. 
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A Equations of State for the Replica-Symmetry-Breaking Solution 

In this appendix we derive the equations of state of the one-step RSB solution 
discussed in section 4. We shall consider the same matrix block sizes, i.e., the 
same value of m for all matrices; denoting the elements of the off-diagonal 

blocks by gi^oj (I2.0: f^o, 7i,0i 72,0 and ,^0i whereas those of the diagonal blocks by 
q'2,1, ''"1, 7i,i, 72,1 and .^1, one gets the free-energy density defined in Eqs. 
(27) - (30). The extremization of such a free-energy density with respect to 
this set of parameters leads to the following equations of state. 
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7i,o=2P(/^Ji)V,o' (A.l) 

72,o = ^p(/5J2)\f,o' (A.3) 

l2,i^\p{(3J2fqlf (A.4) 

io = \p{^J^frl' (A.5) 

^i-\pWrV (A.6) 



?i,o 



J/0, yi, J/2 



'2/0, 2/1, J/2 



Xo,Xl,X2 



'2/0,2/1,2/2 



((A C*)) 2/0, 2/1, J/2 



J/0,yi,2/2 / / a;o,xi,a;2 
n 2\ \ 



pm-2^2 



J/0,2/1,2/2 
2/0,2/1,2/2 



'2/0,2/1,2/2 



a;o,2:i,a;2 



a:o,a;i,a;2 



n 2\ 



J/1,2/2 



2/0,2/1,2/2 
2/0,2/1,2/2 



2/0,2/1,2/2 



a;o,a;i,i;2 



Xo,Xl,X2 



(A.7) 
(A.8) 
(A.9) 
(A.IO) 
(A.11) 
(A.12) 



where we have used the notation introduced in section 3, i.e., the double 
brackets (( ))a;,j/,2 stand for Gaussian averages with respect to the set of variables 
{x, y, z). The extremization with respect to the parameter associated with the 
block sizes, m, leads to 



1 



+ ^ (Ti,i?i,i - 7i,ogi,o + 72,152,1 - 72,o?2,o + 6n - ^oro) 



j3rn/ 



{{H{An)yo,vuv.)) 



(3m 



((A'" In A)) 



Xo,Xl,X2 
1 2\ \ 



2/0,2/1,2/2 



2/1,2/2 



(A.13) 



Xo,Xl,X2 
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In the equations above one has, 



A — A cosh(iio) cosh(iii) cosh(ii2) + 4 sinh(iio) sinh(iii) sinh(ii2) (A. 14) 

B — A sinh(wo) cosh(ui) cosh(ii2) + 4 cosh(iio) sinh(iii) sinh(ii2) (A. 15) 

C = 4 cosh(-uo) sinh(-ui) cosh(-U2) + 4 sinh(-uo) cosh(-ui) smh(-U2) (A. 16) 

D — 4 cosh(Ko) cosh(Hi) sinh(K2) + 4 sinh(Ko) sinh(iii) cosh(K2) (^-l^) 



with 



uo = ^/ylfi xo + 1/71,1 - 7i,o yo 



ui = y/y^ + ^72,1 - 72,0 yi 
U2 



^ X2 + v'6 - 6 y2 ■ 



(A.18) 

(A.19) 
(A.20) 



B Equivalence with Random Energy Model 



In this appendix we consider two particular cases of the model defined in Eq. 
(1), in the limit p ^ 00, and find the corresponding equivalences with random- 
energy models. Let us first consider spin configurations for which {cij} and 
{tj} are completely independent of each other; this is expected to be valid 
at high temperatures, and the corresponding random-energy model should 
yield the same results as the model of Eq. (1) in the paramagnetic phase. The 
probability that such a configuration presents a given energy E is 



P,{E) = {6{E-n{{a„n})), (B.l) 

where the average is taken over all possible realizations of Ji^...i^. Using Eqs. 
(1) and (2) we find 



P,{E) 



exp 



E'' 



(B.2) 



Jl. The 



where we have introduced the effective variance J^^j — Jf 
total number of such configurations is 4^; therefore, the average number of 
macroscopic states with energies between E and E + dE is given by 



ni(£;) = 4^Pi(£;) «exp 



N 21n2 



^2 



(B.3) 
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In the thermodynamic hmit, iV ^ oo, it is convenient to introduce the energy 
u — E/N, and entropy s{u) — S{E)/N densities. From Eq. (B.3) we obtain 



s{u) 



u 



+ 21n2. 



(B.4) 



Using the thermodynamic definition of temperature 1/T — ds/du we may 
write 



s{T) 



"&+21n2. 



4T2 



(B.5) 



When J I = J 2 we recover the entropy density of the paramagnetic phase, 
given by Eq. (22). 

Let us now consider another particular case, namely, situations in which the 
{cTj} and {Tj} configurations are the same; this is expect to hold throughout 
the SGII phase. In this case, the last term in (1) does not contribute, and we 
get the simplified Hamiltonian 



(B.6) 



l<ii<...<ip<N 

The number of such configurations is 2^. The probability that one of these 
configurations presents an energy E is given by 



P2iE) = {6iE~ni{a,})), (B.7) 

with Ti, given by Eq. (B.6). From the above expression and the probability 
distribution for the couplings, we find 



P2{E) 



■■ exp 



^2 



N{J! + Jl) 



Therefore, the mean number of such configurations is 



(B.8) 



n2{E) = 2^P2{E) ^exp{N 



In 2 



(B.9) 



For Jl 

by 



N{J! + Jl) 

J2, the entropy density as a function of the energy density is given 



s{u) = 



u 
"2jf 



+ ln2. 



(B.IO) 



15 



Expressing this entropy as a function of the temperature, one gets the same 
entropy density of the phase SGII [cf. Eq. (35)], 

MT) = -^^+\n2. (B.ll) 

In summary, we have found the equivalence of our model with random-energy 
models in two particular cases: (i) spin configurations {cTj} and {tj} indepen- 
dent from each other. This leads to a random-energy model, appropriated for 
the paramagnetic phase, with 4^ independent random-energy levels Ei follow- 
ing the distribution given by (B.2); (ii) the same spin configurations {cxj} and 
{tj}. In this case one gets a random-energy model, appropriated for the SGII 
phase, with 2^ independent random-energy levels distributed according to 
Eq. (B.8). To obtain the free-energy density for the phase SGII we can use 
the method introduced by Derrida [8] , which avoids the use of replicas. 
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